Nonlinear resonant tunneling in systems coupled to quantum reservoirs 
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An adiabatic approximation in terms of instantaneous resonances is developed to study the 
steady-state and time-dependent transport of interacting electrons in biased resonant tunneling 
heterostructures. The resulting model consists of quantum reservoirs coupled to regions where the 
system is described by nonlinear ordinary differential equations and has a general conceptual interest. 
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The mathematical method recently proposed in ^] 
provides a significant advance in the solution of time- 
dependent scattering problems for Schrodinger equations 
with nonlinearities concentrated near the resonances of 
the corresponding potential. The method consists in the 
separation of the original system in two coupled subsys- 
tems through the splitting of the wavefunction in two 
components. In this way, we can separately study the 
simple problem of a reservoir-like subsystem having only 
extended states and couple its solution to the other sub- 
system having resonance states. The solution of this sec- 
ond problem is then simplified by an adiabatic approxi- 
mation in terms of instantaneous resonances. 

The situation investigated in |Q] depicts ballistic trans- 
port in a double barrier heterostructure, i.e., the scatter- 
ing of a wave-packet on a double barrier potential. The 
nonlinearity, concentrated in the well between the barri- 
ers, is due to the interaction of the electrons represented 
by the wave-packet (mean-field). Here, we generalize the 
approach |l| to the case of biased heterostructures where 
a band of scattering states are to be considered. For these 
systems, hysteresis in the current- voltage characteristics 
has been observed Q and recognized as a consequence 
of the mutual interaction of the electrons trapped in the 
resonance |^],|| . We show that the approach (lj allows us 
to quantitatively reproduce experimental results like Q 
and predict new time-dependent properties. Although 
illustrated in the case of heterostructures, these results 
are very general and have applications in fields like the 
theory of electric systems || and nonlinear optics (6|J^] . 

Let us consider a heterostructure whose conduction 
band edge profile consists of two barriers of height Vq 
located in [a, b] and [c, d] with a < b < c < d along the 
growth direction x. Translational invariance is assumed 
in the plane parallel to the junctions. Suppose that a bias 
energy AV" is applied between the emitter (x < a) and 
collector (x > d) regions uniformly doped with net donor 
concentration np>. At thermal equilibrium with temper- 
ature T, transport is due to a band of scattering states 
with Fermi energy Ep = (3n 2 no) 2 ^ 3 (we use everywhere 
effective atomic units h = 2m* = 1 and e 2 /e — 2a^ 1 , 
where m* is the electron effective mass, e the dielectric 



constant and as = h 2 e/(m*e 2 ) the effective Bohr ra- 
dius). Due to the translational invariance in the plane 
parallel to the junctions, the single-electron scattering 
state at energy E along the x direction is described by 
the one-dimensional Schrodinger equation 

[-id t -0 2 x + V ch (x) + U{4>, x)] <f>(x, t,E)=0 (1) 

where V c \ ) (x) is the step-like conduction band edge profile 
and U(4>, x) takes into account the applied bias and, at 
Hartree level, the electron-electron interaction. Assum- 
ing ideal metallic behavior in the emitter and collector 
regions, i.e., neglecting the formation of accumulation 
and depletion layers, U((j), x) can be obtained as solution 
of the Poisson equation 
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with boundary conditions U(4>,a) — 0, U(<j),d) = — AV. 
The parallel degrees of freedom are considered through 
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where the Heaviside function 0(E) limits the integration 
in (^) to energies above the bottom of the emitter con- 
duction band, E = 0, as correct for Ep <C AV. 

In general, the solution of (§) can not be handled ana- 
lytically. We will suppose that, due to the accumulation 
of electrons in the well with sheet density 

, r(c+d)/2 

s(ct>) - / dEg(E) / dx |<MM,£)| 2 , (4) 

J J(a+b)/2 

ideal metallic behavior in the we'll [b. c] and ideal insulat- 
ing behavior in the barriers [a, b] and [c, d] hold. Then 
the total potential 14b + U in (|l|) is better rewritten as 
V + W where 

V{x) = [V - AV{x - a) /I] l [aM {x) - AV(b - a) /I 
xl[6,c]0) + [Vo - AV{b-a + x-c)/l] 
xl[c,d](aO - AV l [dt+oo[ (x) (5) 

gives the band profile modified by the external bias and 
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W(s,x)=&Tra- B 1 s{(j>) {(x - a)(d - c)/l l [a , b] (x) 
+(b-a)(d-c)/l l [6)C] (s) 
+ (b-a){d-x)/ll [cA {x)} (6) 

depends on the wavefunction </> through the sheet density 
of electrons in the well s (</>). Here £ = b — a + d — c and 
h XOiXl \(x) is 1 if x £ [xq,xi] and otherwise. 

Following Q, we will solve (Q) with the potential (||- 
||) in two steps. Let Vm(x) — V(x) + Vbl[6, c ](^) be the 
potential obtained by filling the well [b, c] . First we solve 

[-idt -d* + V mi {x) + W(s, x)] p(x, t,E)=0 (7) 

and then we look for cf> in the form (j> = p + v where 
u(x, t, E) should solve 

[-id t -dl + V(x) + W(s,x)]i>=V 1[6 >o] (i)A. (8) 

Substituting (Q) with (^]|) corresponds to decomposing 
the original system in quantum reservoirs coherently cou- 
pled to a small subsystem. The wave function p describes 
an electron at energy E which is delocalized in the emit- 
ter and collector regions (reservoirs) and has an expo- 
nentially small probability to be found in the forbidden 
region [a, d]. The wave function v describes the same 
electron in the double barrier region and is driven by the 
value of p in the well [b,c]. Due to the quasi-localization 
of v in [6, c], we have <fi ~ p in the reservoirs and cj> ~ v 
in the well, with an error which is exponentially small in 
the limit of wide barriers JjJ. 

Equation (|7[) can be solved by evaluating the instan- 
taneous eigenstates of the potential Vmi + W. We put 
p(x,t,E) — exp(—iEt)n(x,t,E) and suppose that AV 
and s are slowly varying functions of time so that also 
fj,(x, t, E) is slowly varying in time. In the emitter region 
x < a we take /x(x, t, E) — /x(x, E) as the sum of a left- 
and right-going plane wave at energy E and propagate 
this expression to the adjacent regions by requiring fj, to 
be of class C 1 . For wide barriers we can use a WKB ex- 
pansion for the potential Vmi + W and explicitly evaluate 
fi in the region [6, c] which is of interest for solving (j§). 

Equation (@) can be simplified by developing v into 
the instantaneous eigenstates of the potential V + W and 
keeping only the contributions from the discrete resonant 
states §. In the case of a single resonant state we put 
v(x,t,E) — exp(—iEt)z(t, E)e(s,x) where e(s,x) is the 
(ground) resonant state of the potential V+W with com- 
plex eigenvalue A(s) = Er(s) — iT{s)/2 

[-A(s) - d 2 x + V(x) + W(s, x)] e(s, x) = 0. (9) 

The eigenfunction e(s,x) is of class L 2 on the contour 
7 = (e l9 \ - oo,0] + aj_U [a, d] [j (d + e ie [0, +oo[) for 6» 
conveniently chosen H and satisfies J dx e(s,x) 2 = 1, 
J dx e(s,x) 9 s e(s,x) = 0. Multiplying (^) with e and 
integrating over 7, we get 



d t z(t, E)=i[E- A(s)] z(t, £7) + S(s, £7) (10) 

with the driving term given by B(s, E) = 
iVo f£dx fj,(x,t,E)e(s,x) and the sheet density (Q) re- 
duced, with small error, to 

s = J dEg(E) \z{t,E)\ 2 ^\\z{t)\\ 2 . (11) 

Explicit expressions of A(s) and e(s, x) can be found 
within the same WKB approximation used for evaluat- 
ing fj, [||. For later use we note that Er(s) = E^ + 77s, 
where £7^ = E a - AV(b - a) /I, E a being the (ground) 
eigenstate of the potential Vo [l]-oo,&] (%) + h c ,+oo[( x )\ 
and r] = Sira^fi - a)(d -c)/£= e 2 /(C c + C c ), C c and 
C c being the emitter and collector capacitance per unit 
area, respectively. Moreover, T(s) — T c (s) + T c (s), T c 
and T c being the contributions to the resonance width 
given by the emitter and collector barriers, respectively. 

The original problem (|l|) is reduced to solving the sys- 
tem ( [To| ) with the condition jll]). Let us first consider 
the stationary solutions 

m = v( w^ff p ( w ( 12 ) 

-T\s)/2 + i[E - Er(s)\ 

Equation (|ll]) gives a self-consistency condition for s = 
\\z\\ 2 . Assuming that \B(s,E)\ 2 is a smooth function of 
E and T(s) <C Ep. a Dirac S approximation can be used 
to get 

s = f(s) = 2tt g(E R (s)) \B(s,E R (s))\ 2 T(s)- 1 . (13) 

The function /(s) vanishes everywhere except for 
< Er(s) < Ep where, for Ep <C Vq we have 
2n\B(s,E R (s))\ 2 ~ r c (s). Equation d|) is then equiva- 
lent to sT(s) = g(E R (s))T c (s) which has a simple inter- 
pretation in terms of charge conservation. In the steady- 
state, the current of electrons injected from the emitter 
into the well, g(E R (s))T e (s), equilibrates the escaping 
current, sr(s). The latter current increases with increas- 
ing the sheet density of electrons in the well, s, while the 
former vanishes at both E R (s) — (square root singular- 
ity) and E R (s) ~ Ep. Therefore, Eq. ( p^| ) has only one 
solution for E R > and may have three for E R < 0. For 
E R > Ep the unique solution vanishes and for E R < 
the couple of nonvanishing solutions is to be searched in 
the interval -E R /r) <s< {-E R + E F )/rj. 

In terms of applied bias, multiple solutions of (^) can 
be obtained for > E £/(b - a). The range of AV 
values for which three solutions exist depends on the am- 
plitude of the function /. If the emitter barrier is more 
opaque than the collector one, / is suppressed by the 
factor r c /r c <C 1 and we always have only one solution. 

The solutions of ([l3]) can be characterized in terms of 
stability. This is particularly important in view of a com- 
parison with steady-state experiments where only stable 
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FIG. 1. Theoretical steady-state current-voltage char- 
acteristic for the GaAs-AlGaAs heterostructure experimen- 
tally investigated in under forward (right-most curve) and 
reverse (left-most curve) bias. In the reverse bias case we 
permuted the barriers instead of making AV (and /) nega- 
tive. In the forward bias case the dashed line is an unstable 
steady-state solution and arrows indicate the transition ex- 
pected at the bistability thresholds A and B by decreasing 
or increasing the bias, respectively. The relevant parameters 
are n D = 2 x 10 17 cm" 3 , T=1K, i = 2 x 1(T 5 cm 2 , 
m* = 0.067 m, where m is the free electron mass, e = 11.44, 
V = 0.34 eV, b-a = 9.0 nm, c-b= 5.6 nm, and d-c = 10.7 
nm @. 

solutions are measured. By studying the eigenvalues of 
the linearization of the vector field defined by the r.h.s. 
of ( p~0| ) one can demonstrate that a solution of ( |i~3| ) is 
stable (unstable) when d s f(s) < 1 (> 1) |§. The trivial 
solution s = 0, when it exists, is, therefore, a stable one. 
When three solutions exist, two of them, the largest and 
the smallest one, are stable while the intermediate one is 
unstable. 

Considerations analogous to those made for s hold for 
the steady-state collector current I/e = As T c (s), pro- 
portional to the number of electrons in the well, As, A 
being the transverse area of the heterostructure, and to 
collector escape rate T c (s). In Fig. 1 we show I(AV) 
evaluated for the asymmetric double-barrier heterostruc- 
ture experimentally investigated in In agreement 
with the above discussion and with the experimental find- 
ings, no multiple solutions are obtained in the left-most 
curve (reverse bias case) of Fig. 1 when the emitter bar- 
rier is wider than the collector one. On the other hand, 
a bistability region extending between points A and B is 
observed in the forward bias case. 

The above results can be generalized to include the ef- 
fect of inelastic processes if we change T — > T + Ti in Eq. 
( |l2| ) , where Ti = Ti e + T. kc is the total width representing 
collector and emitter inelastic decay channels [|n]. The 
collector current becomes I/e — As (T c (s) + r,- c ). How- 
ever, if T + Ti < E F , Eq. dT|) still holds so that for 
r c ~ r (Ti c — Ti), as in the case of Fig. 1, the value of 




FIG. 2. Sheet density of electrons in the well s(t) after an 
instantaneous increase SV of the bias from point B of Fig. 1. 

Now we turn to the time-dependent transport proper- 
ties. According to (nffl) and (O), we have 



d t s(t) = -T{s(t)) s(t) + 2Re(z\B) 



(14) 



where (u\v) = J dE g(E) u(E)v(E). The last term in 
( |l4| ) can be expressed in terms of s(t) by using the for- 
mal solution of (|i"c| ) 

z{tjE) ^j:^Hns(t')) +t (E-E R ( S{ t')))] z{0fE)+ 

, e r; d4 ''[-ir( S (0)+^-wo))] B (a(t / )s E) (15) 



dt 



The first term in ( |i~5| ) vanishes exponentially and can be 
neglected after a time t ^> 2/T. In the second term, 
an analogous exponential factor selects the contributions 
for t — t' < 2 /r as the dominant ones so that the lower 
integration bound can be safely changed into -co for 
t 2/T. In this case, 2Re(z\B) can be approximated 
with 



2Re 



dt' Jl d *"[-i r ( s (*"))+^ R ( s (t"))]^ (5 | B |2 )(t _ tl) 



2Re 



dt' e^^^-^T^B^it-t') 



2ng(E R (s(t))) \B (s(t), E R (s(t)))\' 



(16) 



where T(g\B\ 2 ) is the Fourier transform obtained by per- 
forming the energy integral in the scalar product and the 
approximation in the third line is valid for T and dts(t) 
small. With this result Eq. ( |i~4|) becomes 



d t s(t) = -T(s(t)) [s(t) - f(s(t))]- 



(17) 



As Eq. (|T^), this is a conservation law for the charge 
trapped in the well. The steady-state solutions, dts = 0, 
of ( |i"7j ) coincide with those defined by ( |l3| ) and also their 
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FIG. 3. Survival time of the quasi-stationary solutions of 
Fig. 2 versus SV (solid line). The dashed line is Eq. (p^[). 
Arrows indicate the relevant energy scales. 

attractive nature agrees with the stability condition dis- 
cussed above. 

Close to a bistability threshold, the dynamics of cer- 
tain nonlinear optical systems has been shown to be char- 
acterized by a quasi-stationary transient followed by a 
fast evolution As we will discuss in a moment, 

this behavior is typical of any rate equation of the form 
dtS = h(s,SV) where SV is the control parameter of a 
bistability threshold. An example for our heterostruc- 
ture is shown in Fig. 2. At time t = the system is 
in the threshold stable steady-state B of Fig. 1 when 
the bias is instantaneously increased by an amount SV. 
For SV smaller than a critical value, we observe a quasi- 
stationary s(t) which decreases linearly at small times 
and, after a time t, vanishes in a nearly exponential 
way. The border between these two regimes is given 
by the condition E r (s(t)) — 0. Indeed, when the res- 
onant energy E R (s) becomes smaller than the emitter 
band edge the filling current, T(s)f(s), vanishes and jl7| ) 
has solution s(t) cx exp[— f dtT(s(t'))]. The nearly 
exponential decay is established from the beginning if 
SV > E R (s(0)). For SV < E R (a(0)) and t < r, s(t) is in 
a quasi-stationary regime which is characterized only by 
the fact that the starting point, s(0), is a threshold stable 
steady-state. Indeed, for s(t) — s(0) and SV small in this 
case we must have h(s,SV) ~ — q[s(i) — s(0)] 2 — pSV with 
q,p > 0, independently of h(s,SV). Integrating, we get 
s(i) = s(0) — \J SVp/q tan(\/ qpSVt) which has linear be- 
havior for small t. In the case of Eq. ([l7|), the condition 
Er(s(t)) — and the approximate evaluation of q and p 
at T = K give 

where E R (s(0)) ~ E F T C /(2T)[1 + 27rT/(rir e )}~ 1 and 
s(0) ~ [^F-^fl(s(0))]r e /(27rr). When SV < E R (s(0)), 



we have r ~ SV~ X I 2 as shown in Fig. 3 where we com- 
pare ( |l8| ) with r obtained by numerically integrating ( |T^ ) 
|Q. For E R (s(0)) <C Ep, the temperature dependence 
of ( |T8| ) is easily obtained by substituting the Fermi en- 
ergy with the effective value Ep = Ep + fcsTln[l + 
exp(— Ep/k R T)]. The survival time increases by increas- 
ing T and for SV < E R (s(0)) we have r ~ E F /2 . 

The phenomenon discussed above may be exploited for 
device applications like those suggested for optical sys- 
tems 0. More complex time-dependent features are ex- 
pected in heterostructures with many resonances |{J or 
in superlattices 
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